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Abstract. Fix a strictly positive measure W on the d-dimensional torus T'*. 
For an integer > 1, denote by , x = {xi, . . . , xj), < Xi < N, the W- 
measure of the cube [x/N, {x+l)/N), where 1 is the vector with all components 
equal to 1. In dimension 1, we prove that the hydrodynamic behavior of a 
superposition of independent random walks, in which a particle jumps from 
x/N to one of its neighbors at rate {NW^)~^ , is described in the diffusive 
scaling by the linear differential equation dtp = {d/dW){d/dx)p. In dimension 
(i > 1, if IV is a finite discrete measure, W = 5Zi>i '"'i'^i^i i ™s prove that 
the random walk which jumps from x/N uniformly to one of its neighbors at 
rate {W^)~^ has a metastable behavior, as defined in [2], described by the 
/C-process introduced in 13 . 



1. Introduction 

Scaling limits of randoms walks in random trap environments have been ex- 
amined recently [HI El Hj as stochastic models which exhibit aging [HI |3l [7], a 
phenomenon of considerable interest in physics and mathematics. 

To describe the dynamics, fix an unoriented graph G = {V, E) with finite degree 
and consider a sequence of i.i.d. strictly positive random variables {^z : z € V} 
indexed by the vertices. Let {Xt : i > 0} be a continuous time random walk on V 
which waits a mean exponential time at site z, at the end of which it jumps to 
one of its neighbors with uniform probability. 

The time spent by the random walk on a vertex z is proportional to the value of 
^z- It is thus natural to regard the environment as a landscape of valleys or traps 
with depth given by the value of the random variables {£_z ■ z £V}. As the random 
walk evolves, it explores the random landscape, finding deeper and deeper traps, 
and aging appears as a consequence of the longer and longer times the process 
remains at the same vertex. 

It is clear from the description that random walks on random trap environments 
should present a very rich scaling fractal structure if one chooses appropriate graphs 
and random environments. For each given time scale, only traps at a certain depth 
matter. The deeper valleys are too sparse to influence the evolution and the shal- 
lower wells are not deep enough to retain the process. 

We are concerned in this article with the lattice case: {^z '■ z £ U^} is a sequence 
of i.i.d. strictly positive random variables and {Xt : t > 0} a continuous time 
random walk on Z'^ which waits a mean fz exponential time at site z, at the end 
of which it jumps to one of its neighbors with probability l/2d. 

When ^0 has finite mean, for almost all environments : z £ Z''}, the rescaled 
random walk eXt^~2 converges in distribution to a Brownian motion. In dimension 



Key words and phrases. Trap models, scaling limit, hydrodynamic equation, gap diffusions, 
metastability. 



2 



M. JARA, C. LANDIM, A. TEIXEIRA 



1, we can use the method of random time change to study the problem exphcitly and 
a simple computation establishes the result [20j . In this case, the diffusion coefficient 
is equal to E[£^o]~^, the harmonic mean of the random rates {^7^ : z e Z''}. 
Observing that the random walk is a martingale, in higher dimension, by examining 
the evolution of the environment as seen from the position of the random walk, the 
proof of the invariance principle is reduced to the proof of an ergodic theorem 
for the dynamics of the environment [l9j . An explicit formula for the variance is, 
however, no longer available. 

To investigate the case where the environment has an infinite mean, a natural 
assumption is to suppose that the distribution of belongs to the domain of 
attraction of an a-stable law, < a < 1. The variables {i^^ : z E Z'^} take now 
large values in certain sites, forcing the random walk to stay still for a long time 
when it reaches one of them, causing a macroscopic subdiffusive behavior. 

In dimension 1, Pontes, Isopi and Newman [ill proved under these hypotheses 
that for almost all environments, the random walk converges, in the time scale 
^i+(i/q)^ to a singular diffusion with a random discrete speed measure. In dimen- 
sion d > 2, Ben Arous and Cerny ,5,^ proved that for almost all environments the 
Bouchaud trap model converges in a proper time scale, t^^" in dimension c? > 3 
and a scale logarithmic smaller than t^/" in dimension 2, to the fractional-kinetic 
process, a self-similar, non-Markovian, continuous process, obtained as the time 
change of a Brownian motion by the inverse of an independent a-stable subordi- 
nator. In fact, they proved, under quite general conditions on the environment, 
that the clock process converges to an a-stable subordinator, for a large range of 
time scales [6 . In these time scales, the random walk does not visit the deepest 
traps, but exhibit an aging behavior. During the exploration of the random scenery, 
the process discovers deeper and deeper traps which slow down its evolution, the 
mechanism responsible for the aging phenomenon. We refer to ^ i9j for recent 
reviews. 

We present in this article two results. The first one establishes the hydrody- 
namic behavior, almost sure with respect to the environment, of a superposition 
of independent random walks evolving on the one-dimensional torus with a trap 
environment of a-stable i.i.d. random variables. The hydrodynamic equation, de- 
scribing the macroscopic evolution of the density, is given by the generalized second 
order linear equation 

d , d d , 

where W is an a-stable subordinator deriving from the realization of the environ- 
ment. The Krein-Feller operator {d/dW){d/dx) is the generator of the singular 
diffusion obtained by Pontes, Isopi and Newman [llj as scaling limit of the random 
walk in the trap environment. 

The striking feature of this result is that the random environment survives en- 
tirely in the limit, since even the differential operator, which describes the macro- 
scopic evolution of the density, depends on the specific realization of the environ- 
ment. A similar phenomenon was observed in [TUl [T31 [12] for exclusion processes 
with a-stable random conductances. 

The second result describes the evolution of the random walk in the random 
environment, produced by a-stable i.i.d. random variables, in dimension d > 2 in 
the time scale needed to visit the deepest traps. In the notation of Theorem 4.1 in 
[6], this corresponds to the case 7 = 0. 
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In dimension 2, on the time scale iV^/" logiV, we prove that the random walk, 
evolving on the discrete torus (Z/iVZ)^, converges to the Markov if-process intro- 
duced by Pontes and Mathieu [13j . which in the present context can be informally 
described as follows. The state space is formed by the countable and dense subset 
of deepest traps. The process stays at one of these sites an exponential time, with 
expectation proportional to the depth of the trap, at the end of which it jumps 
to a new location, chosen with uniform probability among the deepest traps. The 
scaling limit is similar in dimension d > 3, but the time scale is now iV*/". In the 
terminology of [2j, these results establish the metastability of the random walk in 
dimension d > 2. 

Convergence to the if -process has been proved by Pontes and Mathieu [13j for 
the trap model in the complete graph and by Pontes and Lima [121 for the trap 
model in the hypercube. We believe that this is a universal behavior of random 
walks on graphs with heavy tailed random trap environments in the ergodic time 
scale, the scale proportional to the time needed to jump from one very deep trap 
to another. At least in sufficiently high dimension. 

It is in fact quite surprising that even in low dimensions the geometry of the 
torus is completely wiped out in the scaling limit of the random walk in a random 
trap environment, as proved below. 

We conclude this introduction by specifying the random environment we consider 
in this article. Though we shall work on the torus, we present the construction on 
R'^. Let A be the measure on M"* x (0, oo) given by A = aw~^^^°'^ dx dw, < a < 1. 
Denote by {(x^, Wi) G M"^ x (0, oo) : i > 1} the marks of a Poisson point process of 
intensity A, and define the measure W on by 

i>i 

Por z — (zi, . . . ,Zd) in Z'', let [z/N, [z + 1\/N) be the d-dimensional cube ni<i<(i 
[z^/N, [z^ + l]/N) and let 

i>i 

where 1{A} stands for the indicator of the set A. We show in the next section 
that, for each > 1, {^f^ : z G U^} are i.i.d. random variables with a common 
a-stable distribution, independent of N. Pollowing [TTl Section 3], we may refine 
this construction to obtain i.i.d. random variables distributed according to any law 
in the domain of attraction of an a-stable law. 

Taking the array : z G Z''}, iV > 1, as our environment, instead of a 
sequence {^2 : z G Z"*} of i.i.d. random variables in the domain of attraction of an 
a-stable law, as it is usually done, produces noticeable differences in the scaling 
limit, the main one being the survival of the measure W . 

2. Notation and Results 

Pix a finite, strictly positive measure W on the d-dimensional torus T'': 

W{A) > for any open set A. (2.1) 

Denote by the d-dimensional, discrete torus ['L/N'LY. Let , x G T^, be 
the VF-measure of the d-dimensional cube [x/N, (a; + l)/iV), where 1 is the vector 
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with all components equal tol: 1 = (1,...,1): 

= W{[x/N,{x + 1)IN)] . (2.2) 

We examine in this article the evolution of a continuous time, nearest neighbor, 
symmetric random walk on which waits a mean exponential time at site 
X. Its generator Cn is given by: 

(£iv/)(^) = YdW^ ^t^^^^ ~ -^^"^^^ ' ^^-^^ 

for every / : T^^ R, where {yi,...,yd) = y ~ x = (xi, . . . ,Xd) if \y - x\ = 

J2i<i<d - yi\ = 1- 

2.1. Hydrodynamic limit in dimension 1. Consider a finite number of random 
walks evolving independently on Tjv according to the dynamics defined by the 
generator jCn- Let No be the non- negative integers: No = {0,1,...}. Denote by 
VIn = NJ™ the state space of the process and by 77 the configurations of Clj^ so that 
r]{x), X € Tn, represents the number of particles at site x for the configuration rj. 

This evolution corresponds to a Markov process on fijv whose generator Ln is 
given by 

(Ln/M = ^ E E ^ W") - M ' 

where / : fijv — > R is a bounded function and 77^'^ stands for the configuration 
obtained from r] by moving a particle from site x to site y: 

{r}{x) — 1, z = X 
riiy) + l, z = y 
r]{z), z^x,y. 

Notice that we have slowed down the dynamics by a factor N. We did that 
in order to have a jump rate NW^ of order one if the measure W is absolutely 
continuous with respect to the Lebesgue measure in a neighborhood of x/N. Indeed, 
in this case, if we denote by w the Radon-Nikodym derivative of W, NW^ = 
^ I[x/N {x+i)/N) '^iy)^y is of order one. In contrast, if W has a point mass at x/N, 
NW^ is of order N, which means that particles wait exponential times of order A'' 
at sites where W has point masses. Particles are thus trapped on these sites. 

Denote by {% :t> 0} the Markov process with generator Ljv speeded up by iV^. 
Let I?(R+, r^iv) be the space of right continuous trajectories ^ : M+ fl^ with left 
limits, endowed with the Skorohod topology. For a measure /j. on Qat, let be 
the probability measure on £)(IR+, J]//) induced by the Markov process {rjt ■ t >Q} 
starting from 11. 

For p > 0, let be the Poisson probability distribution with parameter p in 
No: ^p{k} = e~^p'^ /k\, k >Q. Denote by Vp the product measure on 0,^ with 
marginals defined by 

u^{f1 : f]{x) = fc} = «PpM/«{^} , a; € Tat , fc > . (2.4) 

It is not hard to see that the measures I'j^ are invariant and reversible for the 
generator Ljv- 
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Let A^(T) be the space of finite positive measures on the torus T, endowed with 
the weak topology. Fix 7 > and denote by = ir^ [t]) G M{T) the measure 
obtained from a configuration i) by assigning mass N~'^ to each particle: 

'^'^ = 7^ E ^(^)^x/iv, (2.5) 

where ^^/n stands for the Dirac's measure at x/N. For a continuous function 
if : T — > M, denote by (tt^, H) the integral of H with respect to 'k'^ so that 



Fix a continuous function uq : T — > M+ and denote by At^(.) the product measure 
on Qjv with marginals given by 

l^uo{-){'n ■ V{x) = k} = Vuoix/N)Nyw^{k} , X eTN , k>0 . (2.6) 

When uq is constant function equal to p, we denote simply by = /ip. Thus, 

under vi^) ^ Poisson distribution with parameter uq{x/N)W^N^. 

An elementary computation shows that {n^ ,H) converges to / H{x)uo{x)W{dx) 



in L^{iJ,^ , for every continuous function H: 



lim E,,N 



((tt^, i?) - f H{x)uo{x)W{dx)^ ' 



. 



The hydrodynamic equation. Let Tii be the Sobolev space of all functions in 
L^(T) with generalized derivative in L^(T) endowed with the scalar product (•, •)i,2 
defined by 

(/,5)i,2 = {f,g) + / {d^f){x){d^g){x)dx , 
Jt 

where (•, •) stands for the usual scalar product of L^(T). It is well known that the 
space of functions with continuous partial derivatives of all order is dense in Wi. 
Moreover, any function in Hi has a continuous version. 

Denote by L'^{dW) the Hilbert space associated to the measure W{dx), and by 
{f,g)w the corresponding inner product. 

Definition 2.1. A bounded measurable function u : [0,T] x T — > M is a weak 
solution of 

d _l d d 

di^'~2dWdx^' (2.7) 
14(0,-) = uo(-) ; 



if 

(i) It has finite energy: 



I 



T 

{ut,ut)i,2dt < 00 , 
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(ii) For any smooth function G : [0, T] x T — > M vanishing at T, Gt = 0, 
{Go,uo)w + j {dtGt,Ut)wdt = ^ j {dxGt,da;Ut) dt . 
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We prove at the end of this article that there is at most one weak solution of 
()2.7p . Denote by , t > 0, the empirical measure associated to the state of the 
process at time t: 

x6Tiv 

and recall that time has been speeded up by N"^ . 

Theorem 2.2. Let W he a finite, positive measure on T satisfying (|2.ip . Assume 
that there exists 70 > such that 

Fix 7 > 70- Then, for every t > 0, every continuous function iJ : T M, and 
every 5 > Q, 



lim P„iv 



{tt^,H)- / H{x)u{t,x)W{dx) 
where u is the unique weak solution of p.7p . 



> 5 



= 



If the measure W is absolutely continuous with respect to the Lebesgue measure 
and its Radon-Nikodym derivative, denoted by w{x), is strictly positive, w > 
a.s., the previous theorem states that the empirical measure tt^ converges to the 
measure Tr{t,dx) = u(t,x)w{x)dx, whose density u is solution of 



u 

u{0, •) = uo{-) . 

The proof of the hydrodynamic behavior of the empirical measure differs sensibly 
from the usual ones due to the space irregularity of the environment. The lack of 
smoothness is reflected in the dynamics by an erratic time evolution. To overcome 
this issue, we average not only in space but also time, investigating the asymptotic 
behavior of the measure 9Jt^ on [0, T] x T, defined by 







which does not capture space and time discontinuities. 

2.2. Metastable behavior of the trap model in dimension d > 2. Fix a 

finite, strictly positive, atomic measure W on the d-dimensional torus T'': 

i>i 

where {xi : z > 1} is a dense subset of T'* and J2i>i < 

Denote by {wi i > 1} the weights of W in decreasing order so that {wi : i > 
1} = {wi : i > 1} and wi > 11)2 > ■ ■ ■ ■ In case of ties, choose the smallest site 
according to some pre-established order. Let {xi : « > 1} be the position of the 
atoms of W corresponding to the weights {wi : i > 1}: 

i>l i>l 

Recall the definition of given in Denote by {X^^ : t > 0} the ran- 

dom walk on with generator L^. Let _D(M+,T^) be the path space of right 
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continuous trajectories uj : M+ with left limits endowed with the Skorohod 

topology. Denote by , x £ T^, the probability measure on D{M.+ ,T%) induced 
by the Markov process {X^} starting from x. Expectation with respect to is 
denoted by . 

Denote by the unique stationary state of the process {X^ : t > Q}. An 
elementary computation shows that is in fact reversible and given by 

Enumerate according to the weights {VKj^} in decreasing order: 

JIat — \Xi , X2 , ■ ■ ■ , Xj^d } , kV'j.jv ^ I'V'^iv ^ ••• ^ I'V'^iv • 

1 2 

In case of ties, choose the smallest site according to some pre-established order. 
Following [S], we call the sites Xj', j fixed, the very deep traps. These are the 
relevant states of the trap random walk on the scale observed here. 

Since W^(T'') is finite, we may assume that for every M > 0, there exists Nq such 
that Xj e [x^/N - {l/2N)l,x^/N + (1/2A^)1], 1 < j < M, for aU N > Nq. 

To define the trace of the process {X^ : t > 0} on a subset F of T^, let Tj^{t), 
t > 0, F C T%, be the time the process remains in the set F in the interval [0, t]: 

ft 

-F/,\ / If vN 



Jo 



and let S^{t) be the generahzed inverse of {t): 



S^{t) := sup{s > : r^(s) < t} 



It is well known that the process {Xt ' : i > 0} defined by 

is a Markov process with state space F, called the trace of {Xf } on t< . 

Let {Yfc : fc > 0} be the d-dimensional, nearest-neighbor, symmetric, discrete 
time random walk on Z'' starting from the origin. For d > 3, denote by the 
probability that {Y/j} never returns to the origin: 

Vd = Po [Yfc 7^ for all fc > 1 ] . 

Let A^j = {xf , . . . , x^^J, 1 < Af < iV, and denote by {xf : t > 0} the trace 
of the process {X^ : < > 0} on the set A^. The second main result of this article 
states that in dimension d > 3, the trace process {X^''^^} converges, as TV t 00, 
to the random walk on {xi, . . . ,Xm} which waits a mean Wj/vd exponential time 
at Xj and then jumps to {xi : 1 < z < Af } with uniform probability. Note that 
we do not rule out the possibility that the process jumps back to the site where 
it was. To state the result, let : — ^ N, be defined by ^N{xj^) — j and let 
j^N,M ^ ^^(xf ^*^). Clearly, {xf : t > 0} is a Markov process on {1, ... , Af }. 

Theorem 2.3. Fix T > and assume that d> 3. As N 00, the law of {X^'^^ : 
< t < T} converges in distribution to a random walk in {I, ... , M} with generator 
£,M given by 

M 
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Moreover, 



lim limsup max TT;^"''" (T)] 

Af^oo AT^oo 1<J<M ^3 L iv n 



where An,ai = \ ^m- 

In dimension 2 the picture is similar, but the process needs to be speeded up by 
log A^. Denote by {X^ : ^ > 0} the random walk on with generator (log A^)£7v, 
where Cn has been introduced in (|2.3p . Hence, {Xf : t > 0} has the same distri- 
bution as {^/fogiv : i > 0}. 

Denote by P^, x E Tj^, the probability measure on Z?(K+,T^) induced by the 
Markov process {X^ ■ t > 0} starting from x. Expectation with respect to is 
denoted by E^. Denote by {X^ ' : t > 0} the trace of the process {X^ : ^ > 0} 
on the set A^^ and let Xf. ' = ^'jv(Xi ' ). 

Theorem 2.4. Fix T > and assume that d = 2. As AI 1 oo, the law o/{xf^'^^ : 
< t < T} converges in distribution to a random walk in {1, ... , M} with generator 
given by 

Moreover, if we denote by T^"'"^(T) the time spent by the process Xf in the set 
An,m on the time interval [0,i], 

lim limsup max TT^"-" (T)l = 0. 

Af^oo jv-oo i<i<M L iv /J 

We prove in Proposition 16. 191 that in dimension 2 the random walk {X^ : t > 
0} with generator Cn does not leave a very deep trap, staying there indefinitely. 
Therefore, on time scales of order 1 the random walk does not move, and on scales 
of order logiV the geometry is wiped out and the random walk jumps from a very 
deep trap to another one, chosen with uniform probability. 

Recall from jT51 Definition 3.1] the definition of the i^-process, a Markov process 
on N, the one-point compactification of N, characterized by two parameters: c > 
and a sequence {7^ > : z > 1} such that X]i>i7j < While 7^"^ represents 
the rate at which the Markov process leaves i, c is related to the behavior of the 
process at the extra point added in the compactification. 

Denote by {Z^' : t > 0} the Markov process with generator Zm- Fontes and 
Mathieu [13l Lemma 3.11] proved that the process Z^'^ converges, as M | cx), to 
the if-process with parameters c — Q and {wi/vd '■ i > 1}. Next result follows from 
this fact and from Theorem 12.31 



Theorem 2.5. Fix T > and assume that c? ^ 3. There exists a sequence '. 
N > 1}, t 00, such that for any sequence {£n '■ N > 1}, £]y < £*j^, £jv T oO; 
the law of {X/^'^" : < < < T} converges in distribution to the K-process with 
parameters {wi/vd : « > 1} and c — 0. Moreover, 

hm max Ef„ [Tjf (T)] = . 

N^oc l<j<£N ^3 

Of course, a similar statement holds in dimension 2, i.e., for £^ in place of £,m- 
In the terminology of Definition 2.1 in [2], Theorem l2.5l states that in dimension d > 
3 the trap random walk {X[^ : i > 0} is metastable with metastates {xi,X2, ■ . ■} 
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and limit given by the X-process with parameters {lii/vd : i > 1} and c = 0. 
Analogously, in dimension 2, the trap random walk {Xf : t > 0} is metastable 
with metastates {xi,X2, ■ ■ ■} and limit given by the _R'-process with parameters 
{2'Wi/TT : i > 1} and c = 0. 

2.3. Bouchaud's trap model. In this subsection we present an example of a 
random measure W which satisfies almost surely assumption (HI). Fix < a < 1 
and let A be the measure on T'^ x (0, oo) given by A = aw~^^^°'^ dx dw. Since A is a 
positive Radon measure, the Poisson point process T of intensity A is well defined. 
Let {{xi,Wi) : i > 1} he the Poisson marks and define the measure W by 

W = ^ 4i ■ 

i>i 

Note that Vr(T'^) is a.s. finite. On the one hand, the random variable r(T'^ x 
(1, oo)) has finite mean which implies that there are only a finite number of Poisson 
marks on T*^ x [1, oo). On the other hand, X]i>i Wil{wi < 1} has finite expectation. 
Note also that W^(A), W{B) are independent if A and B are disjoints. 

Denote by \A\ the Lebesgue measure of a measurable set A C T^. A simple 
computation shows that the random variable Ty(A) has an a-stable distribution 
for any A with \A\ > 0. In particular, the random measure W is self-similar with 
index a/d in the sense that the distributions of W{PA) and l3'^/"W{A) are the 
same for any /3 e (0, 1) and any measurable set ACT'*. 

We call the random measure W{dx) a d-dimensional subordinator of index a. 
For xGT%, define 

Since W{dx) is self-similar with index a/d, {t^;x G T^} is a sequence of i.i.d. 
random variables with common a-stable distribution — W{T'^) which does not 
depend on N. 

Pontes, Isopi, Newman version of the Bouchaud trap model is the symmetric, 
nearest- neighbor, continuous-time random walk on with generator £jv in which 
is replaced by = N'^^°'W^: 

In dimension 1 , the generator on fi^f corresponding to the superposition of in- 
dependent random walks is given by 

Denote hy {rjj : t > 0} the Markov process with generator L'^ speeded up by 
jyi-l-(i/a) g^j^j denote by Trf the empirical measure associated to the configuration 
r]l by formula (|2.5p . Observe that the time scaling is subdiffusive. 

We show below in (|2.8p that assumption (HI) is in force almost surely. Moreover, 
if the Markov process starts from for some continuous function wq : T — > 

by Theorem 12. 2i for almost all measures W, for all t > 0, the random measure 7rf 
converges in probability to the measure u{t, x)W{dx), where u is the unique weak 
solution of (|2.7p . Note that the noise W survives entirely in the limit, even the 
differential equation depends on W. 
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We conclude this section showing that assumption (HI) is in force for 70 > 
(l/a) — 1. Indeed, with the notation introduced above, assumption (HI) can be 
restated as 

rTTT y ^ — ^ a.s. (2.8) 

It is well known that 1 / has finite moments of any order. Denote by mi the ex- 
pectation of ■ The variance of the previous sum is equal to 7V~f3+27-(2/a)}^2^ 
for some finite constant a^. Therefore, by Chebyshev's inequality, for every e > 0, 



1 ^r- r 1 



Ar{2+7-(i/a)} 



> £ 



^2 



< 



£2jY3+27-(2/a) 



Taking e — N ^ , for S > small enough, it follows from Borel-Cantelli that the 
sum in (|2.8p vanishes a.s. provided 7 > (l/a) — 1. 

3. Proof of the hydrodynamic limit 

In this section we prove Theorem 12.21 Fix T > and denote by A^([0,T] x T) 
the space of finite, positive measures on [0, T] x T, endowed with the weak topology. 
For each > 1, consider the measure 9Jl^ on [0, T] x T defined by 



Hence, if we denote by (( OT^, H)) the integral of a continuous function H : [0, T] x 
T R with respect to OT^, we have that 

Let £'([0, T],Ai) be the space of right continuous trajectories tt : [0, T] — + with 
left limits, endowed with the Skorohod topology. Fix a continuous function uq :T 
M+. Let Sat, iV > 1, be the probability measure on D([0,T],7W) x 7W([0,T] x T) 
induced by the initial distribution A'^(.) the pair ({tt^ : < f < T},97t^): 
Qn = P„« o ({Trf : < i < T}, M^)-^. We prove in Lemma [33] below that the 
sequence {Qn : iV > 1} is tight for the uniform topology in the first variable, and, in 
Subsection l3.3[ that all limit points of the sequence {Qn '■ N >1} are concentrated 
on measures ({ttj : < t < r},9Jt) whose first coordinate is absolutely continuous 
with respect to W , ^{t, dx) = v{t, x)W{dx), and whose density Vt is a weak solution 
of the hydrodynamic equation (|2.7|) . Since, by Theorem 15. H there is at most one 
weak solution, for each < i < T, tt^ converges weakly to v{t,x)W{dx), where vt 
is the unique weak solution of (|2.7p . as claimed in Theorem 12.21 

3.1. Entropy estimates. Recall from [15l Section Al.8] the definition of the rela- 
tive entropy H{X\fi) of a probability measure A with respect to another probability 
measure /i defined on the same space, as well as its explicit formula presented in 
[T51 Theorem Al.8. 3]. An elementary computation shows that there exists a finite 
constant Kq such that 

if(/i.iI(.)lMp) < KoN^ (3.1) 
for all iV > 1. In fact iV^'^i7(/x^^^.j|/Zp) converges to J{uo{x)\og[uo{x)/p] — [uQ{x) — 
p]}W{dx) as N oo. 
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Denote by (•, •)/jp the scalar product of L'^{fip) and denote by the convex and 
lower semicontinuous [HI Corollary Al.10.3] functional defined by 

for all probability densities / with respect to fip (i.e., / > and J fdfip — 1). An 
elementary computation shows that 

In if) = E ' ^^^^^ 

Oi(/) = ^/ {y/fW^) - d^^p . 

By [TSl Theorem Al.9.2], if {S^ ■ t > 0} stands for the semi-group associated to 
the generator N'^Ln, for all t > 0, 

^iv«(.)5f Ia^p) + f I]^{f^)ds < i/^«(.)|A.p) , (3.2) 

Jo 

provided stands for the Radon-Nikodym derivative of j-^S"^ with respect to 
Mp- 

3.2. Attractiveness and coupling estimates. Let {fl,!F,P) be a probability 
space and let ^ be a partial order in fl. We say that a function / : — > R is 
increasing if /(rj) < /(^) whenever r] ^ ^. Let A, fi be two probability measures 
in Q. We say that A is stochastically dominated by /i if J fdX < J fdfi for any 
increasing bounded function / : f2 — > R. An equivalent definition is the following. 
We say that a probability measure A defined in 51 x f2 is a coupling of A and ^ if 
A{A xfl) = \{A), A{fl X A) — ^i{A) for any A £ J^. The measure A is stochastically 
dominated by ^ if there is a coupling A of A and /i such that A((?7, di = 1- 

We say that a stochastic process rjt defined in Q is attractive if for any two 
probability measures Ai ^ A2 there is a process {r/l,!]'^) in x J7 such that rjl is 
distributed as the process rjt with initial distribution /i^ for i — 1,2, and such that 
P{Vt — Vt) ^ 1 fo'' t > 0. We call the process {rj}, f]f) a coupling. 

In fijv, we say that 77 ^ ^ if ri(x) < ^(x) for any x € T^r. For this partial order, 
it is easy to see that T]t is attractive. Indeed, since the state space is finite, it is 
enough to show the existence of a coupling for measures fii concentrated on fixed 
configurations 77% with 77^ ^77^. Define 77^ as the process 774 with initial configuration 
r]^. Then, define the process 774 as a copy of the process Tjt, independent of 77^^, and 
starting from fj, where fj(x) — ifix) — r]^{x). Now define 77^ by taking r]^{x) = 
rj}{x) + fjt{x). Since the motion of different particles is independent, it is clear that 
{ril,rif) is the desired coupling as, by construction, 77^^ ^ 77^ for any t>0. 

In terms of stochastic domination, the definition of attractiveness reads as fol- 
lows. If A^ is stochastically dominated by A^, then Aj is stochastically dominated 
by A( for any time t > 0, where AJ denotes the distribution in JIat of the process rjt 
with initial distribution A*. In particular, we obtain the following inequality, which 
we call the coupling estimate: 

Proposition 3.1. Let , A^ be two probability measure on If X} is stochas- 
tically dominated by \^ , then 

E^i[F(77t)] <E;,2[F(77t)] 
for any t>Q and any bounded increasing Junction i<" : ri^v — s- R. 
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Now we need a criterion to decide whether an initial distribution is stochastically 
dominated by another one. In Mo, consider the canonical ordering. It is easy to 
show that ^Ppj is stochastically dominated by whenever pi < p2- Since the 
measures pp are of product form, pp^ is stochastically dominated by pp^ each time 
Pi < P2- More interesting for us, we have the following. 



. Define 

P = ||wo||oo- Then, stochastically dominated by pp for any N > 0. In 



Proposition 3.2. Fix an initial bounded non-negative profile uq -.T 



particular, 



uo(0 



for any t>0 and for any increasing bounded function F : Ojv — > 

The coupling shows that irf^ , DJl^ converge to measures which are absolutely 
continuous with respect to W, the Lebesgue measures, respectively: 

Lemma 3.3. Every limit point Q* of the sequence Qj\j is concentrated on measures 
Tr{t,dx) = u{t,x)W{dx) (resp. 1Xfl{dt,dx) = v(t, x)dtdx) which are absolutely con- 
tinuous with respect to W (resp. the Lebesgue measure) and whose density u{t,x) 
(resp. v{t,x)) is positive and bounded by \\uo\\oo- 

Proof. Fix a limit point Q* of the sequence Qn and assume, without loss of general- 
ity, that Qn converges to Q* (in the uniform topology on the first coordinate). Fix 
a continuous, positive function G : [0, T] x T ^ M, e > and recall that p = HuoHoo- 
By the previous proposition. 



[m,G] 



> 



dt / G{t, x)dx + £ 



< 



■ lip 



T 
T 



;S!Jt,G)) > p f dt f G{t,x)d. 
Jo Jt 



' x -\- e 



for every > 1. We may replace the integral JjG{t,x)du by the Riemann sum 
because G is continuous. Thus, for N large enough, the previous expression is 
bounded above by 

G{t,x/N){^~pN''}dt > e/2 . 



xeTis 



By Chebyshcv and by Schwarz inequalities, since pp is a stationary state given by 
a product of Poisson measures, this expression is less than or equal to 



4T 

72" 



Jo 



dt . 



In view of assumption (HI), this expression vanishes as t oo because G is a 
continuous bounded function. 

Since Qjv converges to Q*, for every e > 0, 

f-T 



Q* 



[^,G)) > p ( dt ( G{t,x)dx + t 

Jo Jf 



. 



Letting e | 0, we conclude that Q* is concentrated on measures 371 such that 
((9Jl, G)) < p Jq dt J.j.G{t,x) dx. Taking a set {Gk : fc > 1} of positive, bounded. 
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continuous functions dense for the uniform topology, we conclude that Q* is concen- 
trated on absolutely continuous measures dx) — v{t, x)dt dx, whose density 
v{t, x) is bounded by p. 

A similar coupling argument shows that for every < t < T and every continu- 
ous, positive function H : T ^ M., 



lim P„2v 



(7rf,i7) > p / H{x)W{dx)+s 







Since we assumed compactness in the uniform topology, we deduce from this formula 
that 



{■^t,H) > p H{x)W{dx)+e 







It remains to recall the arguments presented for DJl to conclude the proof. 
3.3. Hydrodynamic limit. We prove in this subsection Theorem 12.21 



□ 



Theorem 3.4. The sequence of probability measures {Qn '■ N > 1} converges to the 
measure Q* concentrated on the absolutely continuous pair {{nt : < t < r},9Jt), 
TTt — v{t,x)W{dx), 071 = v{t,x) dtdx, whose density v{t,x) is the weak solution of 
the equation (|2.7p . 



Proof. By Lemma [3751 below, the sequence {Qn ■ N > 1} is tight. Fix a limit point 
Q* and assume, without loss of generality, that Qn converges to Q* . 

Fix a smooth function : [0,T] x T ^ R such that H{T, •) = 0. Consider the 
martingale M^{t) defined by 



Mf(<) = {7r^,Ht) - {7r^,Ho) - / {n^,d,H,)ds - / LN{7r^,H)ds 



(3.3) 



The variance of this martingale is equal to 



N 
2N^f 



E E 

x£Tm y:\y~x\ = l 



Vsjx) 



{H{s,y/N)-H{s,x/N)yds 



The coupling estimate shows that the expectation of this expression with respect 



to P„ 



is bounded by CqA^ for some finite constant Co which depends on H 



and p. On the other hand, an elementary computation shows that 



Ln{tt^,H) ds 



1 



m^,ANH) 



where An stands for the discrete Laplacian. In particular, in view of (|3.3p and 
since H{T, •) vanishes, for every S > 0, 

fT 



lim P„iv 



n':,dsH,)ds + {l/2){{m'\ANH 



> 6 







The first term of this sum converges to jj Ho{x)uo{x)W (dx) in P^n -probabil 



o(-) 



ity, as iV I oo. The last expression can be written, up to smaller order terms, as 
(l/2)((9Jl^, A/f )). Hence, since converges to Q*, for every S > 0, and every 
smooth function H , 



Hoix)uo{x)Widx) + / {Trs,dsH,)ds + {l/2){{m, AH)) 



> S 
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Letting (5 J, 0, by Lemma [3T3l Q* almost surely, 

Ha{x)uo{x)W{dx) + [ ds [ {d,H){s,x)u{s,x)W{dx) 



+ (1/2) / ds iAH){s,x)v{s,x)dx = 



According to Lemma l4.61 we may replace u by w in the second term. By Proposition 
14. 3[ we may integrate by parts the last term to obtain that 



{Ho,uo)w + / {dsHs,Us)w ds ~~ (1/2) / {dxHs,dxVs 

Jo JO 



ds = 



This proves that Q* is concentrated on weak solutions of (|2.7p . By Proposition 
4.3i dxV belongs to L'^{[0,T] x T) and by Lemma [3.31 v is positive and bounded. 
Since the previous identity holds for all smooth functions H, v is a weak solution 
of (123. □ 



Theorem 12.21 follows from this result and the tightness in the uniform topology 
of the sequence {Qn ■ N > 1} proved in Lemma [33] below. 

Lemma 3.5. The sequence {Qn '■ N > 1} is tight in the uniform topology in the 
first coordinate. 

Proof. To prove tightness of the sequence {Qn : N > 1} we need to examine the 
two coordinates separately. 

Clearly, the sequence of random measures dJl^ is tight if and only if the sequence 
of random variables (( Tl^ , G )) is tight for every continuous function G : [0, T] x T — > 
M. Tightness of the sequence {{Tl^ , G)) follows from a coupling argument similar 
to the one used in the proof of Lemma 13.31 

To prove tightness of the sequence of processes {tt^ : < t < T} in the uni- 
form topology, it is enough to examine the process (tt^, H) for some fixed smooth 



function H. Recall the definition of the martingale M^{t) introduced in (|3.3p 
Tightness of {tt^, H) follows from tightness of the martingale M^{t) and tightness 
of the additive functional /J N^L]y{'K^ , H) ds. 

The martingale is tight in the uniform topology because, by Doob inequality and 
by the explicit computation of the quadratic variation of M^{t), for every S > 



lim 



sup 

0<t<T 



\M§{t)\ > S 







On the other hand, computing N'^L^{tt^ , H), by Chebyshev and Schwarz inequal- 
ities, for every 5 > 0, 

i-t 

N'^LN{TT^,H)dr\ > S 



sup 

0<|t-s|<is 



< 



eCo 
6^ 



"o(-) 



we may replace the measure by the stationary measure fi^ , estimating the 



for some finite constant Cq which depends only on H . By the coupling estimate, 
we may replace 
expectation by 
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By assumption (HI), this expression is bounded uniformly in N, which concludes 
the proof of tightness. □ 

4. Entropy estimates 

We prove in this section the main estimates needed in the proof of hydrodynamic 
limit. For i> 1, let be a cube of length £: Ai — {1, . . . ,£} and let A^/ = a; + . 
Denote by M^{x) the number of particles on A^j and by Wn{x, t) the W-measure 
of the cube A^/ rescaled by N: 

yGAe y&Ae 

Note that WNix^eN) - W,{x). 

4.1. Two blocks estimate. We prove in this subsection the so called two blocks 
estimate. 



Lemma 4.1. Fix a bounded function G : [0,T] x T ~. 

' 'Vs{x) 



lim lim sup E w 



Mf{x) 



WN{x,tN) 



ds 







Proof. Fix a bounded function G : [0, T] x T ^ M, 5 > and a positive constant 
Ci — Gi{S) to be specified later. Let 



WNix,eN)} ' 



Re{s,vi) = ^ G{s,x/N) ^■ 



By the coupling estimate, 



Re{s,7Js)ds 



y=0 "x+y 



< Cqc^T 



(4.1) 



for some finite constant Cq depending only on Ci, G, p. It is therefore enough to 
prove that 



lim lim sup 



"o(-) 



Reis,Vs)ds 



= 



By the entropy inequality, Jensen inequality and the entropy estimate p.ip , the 
previous expectation is bounded above by 



A 



1 



AN-r 



log E 



expATVjl f V,°{s,r,s)ds 



R. 



(s, 77s) dsj 



for every positive ^ > 0. 

Let A — KqS^^ . Since e'^l < + e^^ and since lim sup„„,ojj iV"''' log {a^y + a^} 
— maxi=i_2 limsup„^oQ N^'^ logajy, to prove the lemma it is enough to show that 
for every S > 0, 



lim lim sup —ttt- log 



e^^^AN^ V;(s,77,)ds} 



< 0, 



where V, = K° - R^ 
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By classical arguments, relying on Feynman-Kac's formula (cf. [El p. 267]), the 
previous expectation is bounded above by 

sup I / V,{s,i^)f{i^)iip{dii)-^I^(J)]ds, 



/ 



where the supremum is carried over all densities / with respect to ^p. Hence, to 
conclude the proof of the lemma, it is enough to show that 

V^{s,ri)f{r,)^,p{d^) < I R,is,rj)f{r^)f,p{dv) + -^I^if) (4.2) 

for every density function / and every S > 0. 

Recall the definition of Wn{x, eN) to rewrite V^{s, rj) as 

Fix a density / with respect to /ip. Performing a simple change of variables, we see 
that 

r]{x) r]{x + y)\ s--^ f r f]{x + z) r]{x + z + 1) 



Z—Q'^ X-\-Z **X-\~Z-\-l 

z=0'' ^^^+^ 



y-1 



Since (a — b) ~ (-^0 — V5)(\/a + Vfo), by Schwarz inequality, the previous expression 
is less than or equal to 



^ 2=0 2 = 0"^ ^ + ^ 

for all /3 > 0. The same change of variables permit to estimate the second term as 



It follows from the previous estimates that for any density / with respect to /ip, 
and aU (3 > 0, 



j K\s,,)m,Mv) < ^ E E w^E/--^ 



+i(/) (4.3) 



We examine each term on the right hand side separately. Set (3 = 2eN~^A. 
Changing the order of summation, we obtain that the second term is less than or 
equal to 

AeA ^ 2 f^vix + y) 



sSTjv y=o " '■^+v 



QUENCHED SCALING LIMITS OF TRAP MODELS 



17 



This expression is bounded by the expectation of with respect to f{ri) d^p pro- 
vide we choose Ci > = AKqS^"^ . By similar reasons, the first term on the right 
hand side of (|4.3p is bounded above by 



iV(eiV+ 1) 
2Aem 



Hence, (|4.2p holds and the lemma is proved. 



□ 



Consider a sequence {Ga^.e : -/V > 1, e > 0} of functions GN,<i ■ [0, T] x Tjv K- 
In the proof of the two blocks estimate, the boundedness assumption on G was used 
only at (14. ip . In particular, the proof presented above shows that 



lim lim sup 



E„ 



"o(-) 



provided 



rT ^2 

limlimsup / ^7 GNe{s,x)'^ds = 



(4.4) 



Recall that : T ^ M is defined by W^{x) = W{[x, x + e]). 
Corollary 4.2. Let J : [0,T] x T R be a continuous function. Then, 



lim lim sup E 



Proof. Since J is a continuous function. 



= 



-i/eA_ 



is absolutely bounded by C{e)N^^ X^xeTjv ^('^) '^'-'^ some finite constant C(e) which 
vanishes as e | 0. In particular, by the usual coupling estimate and changing the 
order of summation, we get that 



lim sup 



JV>1 



"o(-) 



T ^-1 



J{s,x/N)Mf{x)ds 



= 



The second term inside the absolute value can be rewritten as 

Mfjx) 
WN{x,eN) 



l:^ ^ Jis,x/N)e-'WNix,eN) 



7V1+7 



Let GN,eis,x/N) — J{s,x/N)e ^WN{x,eN). Since J is a bounded function, by 
definition of Wn{x, eN), 

E G^A^^x^ds < ^ E wM^r 
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for some finite constant Co which depends only on J. As t oo this expression 
converges to 

CoT [ W2e{xfdx. 



lim lim sup 



For Lebesgue almost all x, W2t{x) vanishes as e | 0. Therefore, by (|4.4I 

Mf{x) ri,{x) 



E 



T ^-1 



Ai+^ 



J{s,x/N)WN{x,eN)[ 



WN{x,eN) 



ds 



Finally, since 



^ \wM{x,eN)-W,{x/N) 



vanishes as A t oo, we may replace Wm{x, eA) by We{x/N) to conclude the proof 
of the corollary. □ 



4.2. Energy estimate. We proved in Lemma [3.31 that any limit point Q* of the 
sequence {QAr:A>l}is concentrated on measures 9Jl which are absolutely 
continuous with respect to the Lebesgue measures: 971 — v{t,x)dtdx. We show in 
this section that the density v has a generalized space derivative in i^([0, T] x T). 

Proposition 4.3. Any limit point Q* oj the sequence {Qn : A > 1} is concentrated 
on measures 971 — v{t, x)dt dx with the property that there exists a function F in 
i2([0,r] X T) such that 



ds / {dxH){s,x)v{s,x) dx = ds H{s,x)F{s,x) dx 

Jt Jo Ji 

for all smooth functions H. We denote the generalized derivative F of v by dxV. 

The proof of this proposition relies on the following estimate. 
Lemma 4.4. Fix a set of smooth functions : [0, T] x T ^ R, 1 < i < i. Let 

' rj{x) r]{x + 1) ~l 



V^is,v) 



Then, for any /? > 0, 



^ E H.is,x/N){ 



lim sup E n 



max 
i<i<; 



2/3 



j:H.is,x/Nr!I^}ds 



< 



Ko 



where Kq is the constant given by (|3.ip . 
Proof. Fix /3 > and let 

A summation by parts and a coupling estimate similar to the one used in the proof 
of Lemma 13.31 shows that it is enough to prove that 



lim sup E 



N- 



max / Xj(s,ri) ds 
i<i<i ' ^ ' " 



< 



Ko 
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By the entropy inequality, Jensen inequality and the entropy estimate p.ip . 

i-T 



"0(0 



< 



max / Xi(s,ris)ds 
i<i<e ' 







1 



Ko 

(3 ^ (3^ 



log 



exp • 



I max PN'^ [ Xi{s,'qs)ds\ 



for every /3 > 0. 

Since, on the one hand, exp{maxi<i<£ a]y} < X]i<j<£ 6xp{a]Y} and, on the other 
hand, limsup„^^ A^"'' log {Y.i<i<ib^N] = maxi<i<£ limsup„_^oo A^"'' logb^r, to 
prove the lemma it is enough to show that 

fT 



limsup-^log E^Jexpl/JTV / Xi{s,?]s)ds\ 



< 



(4.5) 



for 1 < i < £ and any f3 > 0. 

By classical arguments, relying on Feynman-Kac's formula (cf. [El p. 267]), the 
previous expectation is bounded above by 



sup 
/ 



Xi{s,ri)f{ri)^j.p{dr]) 



f3Nf 



tW 



(/)} 



ds 



where the supremum is carried over all densities / with respect to /ip. 
Therefore, to conclude the proof of the lemma, it is enough to show that 



r ri{x) rj{x + 1) 



for all density / and f3 > 0. 

Recall the definition of Vi. Performing a simple change of variables, we see that 



ri{x) ri{x + l) 



{/(^)"/K'^+S)}dA*p 



Since (a — b) — {^/a — \/b){^/a + Vb), by Schwarz inequality, the previous expression 
is less than or equal to 

A f r]{x) 



d^x + lif) 



J ^{^/fi^^^)+^)ydn. 



for all A > 0. The same change of variables permit to estimate the second term as 



A 



r]{x) ri{x + l) 



} fin) diJ-p ■ 



Choosing A = PN^^\H{s,x/N)\, we obtain that for any density / with respect 
to ^p. 



Vi{s,r)){s,ri)f{'r])fip{d'q) < 



which proves the lemma. 



ri{x) rj{x + 1) 



} fiv) dfJ-p , 



□ 
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Recall that, by Lemma l373t Q* is concentrated on absolutely continuous measures 
m = v{t,x)dtdx. 

Corollary 4.5. Any limit point Q* of the sequence {Qn : N > 1} is concentrated 
on measures 9Jl — v{t,x)dtdx such that 

i-T 



sup 

H 



ds / {dxH){s, x)v{s, x) dx 



ds I H{s,x 

T 



Yv{s, x) dx^ 



In this formula the supremum is taken over all functions H in C ■ ([0, T] x T). 

Proof. Fix a limit point Q* of the sequence Qn and assume, without loss of gen- 
erality, that Qn converges to Q* . Consider a sequence {Hj : j > 1} of functions in 
C°'^([0, T] X T) dense for the uniform topology. It follows from Lemma [44l with (3 — 
1, a summation by parts and a coupling estimate, similar to the one used in the proof 
of Lemma [3731 to replace the discrete derivative N{H{s, {x + l)/N) — H{s, x/N)} 
by the continuous one {dxH){s,x/N), that 



max 

l<i<l 



ds / {dxHi){s,x)v{s,x)dx 



ds 



Hi{s, xY v{s, x) o?x| 



< Ko 



Letting £ '\ oo, we conclude the proof of the lemma applying the monotone conver- 
gence theorem. □ 

Proof of Proposition \4.3\ The proof is similar to the one of 15, Theorem 5.7.1] and 
left to the reader. Note that we have in fact 



ds 



{dor:V)is,x)'^ 

v{s, x) 



dx < oo . 



□ 



4.3. Tit ~ T^t, W almost surely. We prove in this section that 9Jlf = nt, W 
almost surely. 

Lemma 4.6. Every limit point Q* of the sequence Qn is concentrated on measures 
Tl{dt, dx) = v(t, x)dtdx and 'K{t, dx) = u{t, x)W{dx) such that u = v (dt x W{dx)) 
almost surely on [0,T] x T. 

Proof. Fix a limit point Q* of the sequence Qn and assume, without loss of gener- 
ality, that Qn converges to Q* . Fix a continuous function J : [0,r] x T ^ M. By 
Corollary [42] and by Lemma [3731 

T 

ds 



J(s, x) u(s, x) W{dx) 

■ T 



lim Eq> 



It follows from the energy estimate stated in Proposition 14.31 that 



ds I J{s,x)e ^ We{x) v{s, x) dx 



lim 



ds I dx J(s, x) — 



{v{s,x) -v{s,y)}W{dy) = 



x,x-\-e) 
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Q* almost surely. Changing the order of summation, it follows from the continuity 
of J that 

lim f ds f dxJ{s,x)- f v{s,y)W{dy) 

•s^OJo Jt e J[x,x+e) 

ds / J{s,x)v{s,x)W{dx) . 



Hence, for all continuous function J, Q* almost surely 

ds / J{s,x) {u{s,x) — v{s,x)}W{dx) — 0, 

JT 

which proves the lemma. □ 

5. Uniqueness of weak solutions 

Theorem 5.1. There exists at most one weak solution of \2. ?[ ). 

Proof. We use a method due to Oleinik (cf. pg. 90 in [H]). Due to the linearity 
of problem (|2.7p . it is enough to show that the constant function equal to is the 
unique weak solution of equation (j2.7p with initial condition uq = 0. 

Fix such solution u. By condition (i), u belongs to L^([0, T]; Tii). Since u{t, •) is 
continuous for almost all t, it is not difficult to show that there exists a sequence 
of smooth functions : [0, T] x T ^ R, e > 0, such that | lie \\00 ^ CXD and 

lim / dt\\\u,{t,-)~u{t,-)\\lw + \\{d^u,)[t,-) - {d^u){t,-)\\l^ - 0. 



C-+0 







Consider the test function Ge : [0, T] x T M defined by 

r 

Ge{t,x) = u^{s,x)ds. 

Since dtGe = u^, and since u^{t, •) converges to u{t, •) in L'^{dW) for almost all t, 
by the dominated convergence theorem, 

lim / {dtGe{t,-),u{t,-))w dt = / {ut,Ut)wdt. 
'^°Jo Jo 

On the other hand, since 

' {{d,G,){t,-),{d,u){t,-))dt = - [ dt [ {{d,u,){s,-)Ad.u){t,-))ds , 



and since S^We converges to dxU in L^([0,T] x T), 



Ihn^ {id,G,)it,-),idxu){t,-))dt = - ^ J^(^J^ {dxu){t,x)dty dx . 
Hence, by condition (ii), since uq = 0, 

J {ut,ut)wdt = J {J {dxu){t,x)dt^ dx . 

This show that ut = for almost every t, and uniqueness follows. □ 

6. Atomic trap models in dimension d>2 
We prove in this section Theorems 12.31 12.41 and [ 



, 2 
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6.1. Capacity and trace process. To help the reader to foUow the arguments 
of this section, we summarize below known results on capacity and trace processes 
used later. Consider a reversible, ergodic Markov chain {Xt : t > 0} on a countable 
set E. Fix a non-empty subset F oi E and denote by {X^ : t > 0} the trace process 
of {Xt : i > 0} on F, as defined in Subsection 12.21 

Denote by v the unique invariant probability measure of {Xt : t > 0} and by 
the invariant probability measure of the trace process {X^ : t > 0}. By Lemma 
5.3 of ^2j, coincides with the measure v conditioned to and is reversible. 

For X £ E (resp. x S F), let P^, (resp. ) be the distribution on the path 
space £'(K+,i?) (resp. D(R+,F)) induced by the process {Xt : i > 0} (resp. 
{Xt : t > 0}) starting from x. 

For a subset B oi E (or F), denote by H{B) the entry time in B, defined as 

H{B) = inf{t > : Zi e B}, 

where Zt stands either for Xt or for X[ . The context will always clarify to which 
process we are referring to. Denote by t{B) the time of first return of {Xt : t > 0} 
to B: 

t{B) = inf{t >Ti:Xte B}, 
where Ti stands for the time of the first jump of {Xt : t > 0}. When the set B is 
a singleton {x}, we denote H{{x}), t{{x}) by H{x), t{x), respectively. 

Denote hy X : E ^ K+ the holding times of the Markov process {Xt t > 0}. 
By Lemma 5.4 in [2], the rate r^{x,y) at which the trace process {Xf : t > 0} 
jumps from a site a; G F to a site y (z F, y x, is given by 

r^{x, y) = \{x) P, [H{y) < t{F \ {y})] . (6.1) 

The expectation of an entry time has a simple expression in terms of the ca- 
pacities associated to the process {Xt : t > 0}. Denote by L the generator of the 
process {Xt : t > 0}. Let A, B C E he two disjoint sets. Define 

B{A, B) ^ {f : E R : f{x) = 1 for a; e A and f{x) = for x e B}. 

Let D be the Dirichlet form associated to {Xt : t > 0}: D{f) = — j JLfdv for 
any f : E ~^ M.. The capacity cap(A, B) between A and B is defined as 

cap(A, B) = ini{D{f) : / e B{A, B)} . (6.2) 

Notice that cap(A, i?) = cap{B,A); it is enough to consider f — 1 — /. An 
elementary computation shows that cap(A, B) = D(fA,B), where /a,_b is the unique 
solution of 

(L/)(x) =0 if a; e F \ (AUB) 
f{x) = 1 ii X e A 

f{x)^0 ifxeB. 
It is easy to see, by the strong Markov property, that 

fAAx)=^AH{A) <HiB)]. (6.3) 

Let A, B be two disjoint subsets of F. Define ca.pp{A, B) as the capacity between 
A and B, with respect to the trace process {Xf : t > 0}. By Lemma 5.4 (d) in [2], 

ce^MAB) = (6.4) 

The first result of this section establishes the relation between capacity and 
expectation of hitting times. 
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Lemma 6.1. For any subset A of F and any y in F \A, 

Proof. Fix Ad F and y \n F\A. By equation (4.12) in [2J, 

K[HiA)] ^ ^£^^Y./i^)nAHiy)<HiA)]. 

To prove the lemma, it remains to recall that is the measure v conditioned to 
F, ((63) and Lemma 5.4 (a) in [2]. □ 

Note that Tz[H{y) < H{A)] vanishes for z in A. In particular, 

Capacities are also related to return times. Next result follows from equation 
(6.10) in [2]. 

Lemma 6.2. Let A be a finite subset of E, and let y G E \ A. Then, 

cap{A, {y}) 



Fy[H{A)<T{y)] = 



When the Markov chain {Xt : i > 0} is not crgodic, the definition of the capacity 
can be generalized in a natural way. Assume that there exists a positive measure 
V, reversible and invariant for {Xt : < > 0}. In this case, of course, v{E) = +oo. 

To define the capacity between a finite set A and infinity, consider an increasing 
sequence of finite sets ^ E such that U„i3„ = E. Since A is finite, cap(^, B^J, 
given by the variational formula (|6.2p . is well defined for n large enough. The 
sequence of functions fA,B^, introduced in (|6.3p . is increasing and bounded. There- 
fore, we can define /yi(a;) = lim„ //i^sc (x). It is not difficult to check that 

fA{x) = P.AH{A) < oo] , D{fA) = inf{D{f) : / e B{A)} , 

where B{A) is the set of finitely supported functions / : _E ^ R such that f{x) = 1 
for X E A. Let cap(A) := D{fA) be the capacity of A with respect to infinity. 

By the dominated convergence theorem and Lemma 16. 2[ the following result 
holds. 

Lemma 6.3. For any y in E, 

cap(y) 



We conclude this subsection with two estimates for the simple symmetric random 
walk on the torus or on the lattice Z'^. Taking adavantage of the commuting 
time identity, which relates expectation of hitting times with capacities. Proposition 
10.13 of [17] establishes the following bounds for the hitting times of the simple 
symmetric random walk on T^: 

Lemma 6.4. Let x, y two points at distance k on the torus T^. There exist 
constants < Cd < Cd < +oo such that in dimension d > 3, 

CdN'^ < E4H{y)] < CdN'^ uniformly in k 
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and in dimension 2, 

C2N^\ogk < E,[H{y)] < C2NHogik + l). 

For > 0, let us denote by the cube of length 2N + 1 centered at the origin: 
— {—N, . . . , N}'^. Denote by SAn its inner boundary: SAn — {x e A^ : 3y ^ 
A^ with \y — x\ = 1}. The following lemma is Proposition 2.2.2 of |16j : 

Lemma 6.5. Let X{t) be the simple symmetric random walk in , d > 3. Then, 
there exist constants < Cd < C'd < +oo such that for any set A C A^ and any 
X € c^AjjY have 

CdN'^^'^cMA) < F^[Ha<+oo] < Cd N^^"^ cap{A) . 

Note that cap(A) is finite because the random walk is transient. 

6.2. Random walks in T^, d > 3. In this subsection we prove some properties 
of the simple random walk on which will be used to establish its metastable 
behavior. Denote by {Yi^ : k > 0} the discrete time, nearest-neighbor, symmetric 
random walk on and let Q^, x e T^, be the measure on I?(Z+, T^) induced by 
the random walk Y'^ starting from x. Expectation with respect to is denoted 
by the same symbol. 

For a subset B of T^, denote by H{B) the entry time in _B, defined as 

H{B) = inf{fc > : e 5} . 

In Lemma 16.61 below we prove that whenever the simple random walk starts 
from a point isolated from the very deep traps, asymptotically, the next very deep 
trap to be visited is uniformly chosen. In Corollary 16.71 we obtain the limiting 
distribution of the next very deep trap to be visited starting from another very 
deep trap. CoroUarv 16.81 presents the limit of the capacity between two points of 
far apart. 

Let d{x,y) be the distance induced by the graph T^. For a subset F of T'' and 
r > 0, denote by B{r,r) the set of sites in at distance less than or equal to r 
from F: B{r,r) = {x G : d{x,r) < r}; and denote by dG the sites not in F 
which are at distance one from F: dT — {x ^ T : d{x,T) = 1}. In these definitions 
we identified F with its immersion in T^: F = F n T^. 



Lemma 6.6. Suppose Ij^ ] oo and A^^ = {xf^ , . . . , xjj;^} C are such that, if 
i j, d{xf ,x^) > In. Then, defining A^^ ^ to be A^^ \ {xf}, we have 

1 



lim sup 



}^[Hix^)<H{AZ,,)]^^ 







Proof. The proof is based on the fact that a site is reached on the scale A'"'^ and 
equilibrium is reached on the scale iV^. Hence, in an intermediate scale, the process 
has not reached A^ and is in equilibrium. In particular, it has a probability 1/M 
to attain x^ before the set A^j 

First, we prove that the process does not reach a site in the scale N^^^: 

lim sup Q^[H{0) < N^/^] =0. (6.6) 



y-d{y,0)>lN 
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By the strong Markov property, 

Q^[H{0) < iV5/2] < [i/(0) < HidBiO,N/8))] 
+ sup 

z<£dB{Q,N/a) 

Since d{y,0) > In, by Lemma 16.51 the first term in the sum above is bounded by 
C'o/^aT^ for some finite constant Cq independent of N. We can therefore suppose 
in dESl) that In = N/8. 

Denote by {Rk ■ k > 1} and {Dk : fc > 1} the successive return and departure 
times between Bi = B{0, N/8) and B2 = B{0, 7V/4): 

Ri = Di = Ri + Hb^^ o ^j?! 

Rn = Dn-i + He, o eD„_, D„ = Rn + Hg^^ oOr^ , n>2. 

Here 0s : D{1+,T%) ^ D{Z+,T%) is the shift map given by YtoO^^ Ys+t- 
For y such that d{y, 0) > by the strong Markov property, 

-9/4 < N^/^] + Q^[H{0) < Rn^-9,.] 

5f [i?^r.-9/4 < N^/^] 



l^[H{0)<N'/^] < Q^[Rn^- 



< sup IS 

z-d(z,0)>N/S 

+ iV^-^/* sup ( 

z£dB{0,N/8) 



(6.7) 



[H{0) < H{dB{0,N/A))] 



The right hand side does not depend on the choice of y and tends to zero a.s N 1 00, 
by [Hj, Proposition 1.1 with u — N~^/^, and Lemma [631 This proves (16.61) . 

In the time scale N^^^ the process reaches equilibrium. More precisely, denote 
by TT^ the uniform probabihty measure on and by — the total variation 
distance between two measures /i, v on T^. By Corollary 5.3 and equation (5.9) 
in [17], for an arbitrary sequence y'^ £ B{A^j,InY , 



lim 



yN \Yn--/ 



(6.8) 



In particular, for 1 < i < M and for an arbitrary sequence y^ in B{A^j^ InY, 



= 



(6.9) 



To prove this claim, fix 1 < i < Af and introduce the indicators of the sets H{xf) < 
iV5/2^ H{xf) > N^/^, to obtain that 



R 



N 



where the remainder Rn is absolutely bounded by 



sup 

zGTi 



Qf[7f(xf)]Qf„[i/(xf)<7V5/2] 



Hence, 

is absolutely bounded by 



R 



N 



N'/' + sup Qf[H{x^)]{Q^4H{x^) < N^/'] 
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By (|6.6p . (|6.8|) and Lemma [6^ this expression divided by N'^ vanishes as t oo. 
This proves ([Q]) . 

Recall that Vd stands for the probability that a symmetric, nearest-neighbor 
random walk on Z"^ never returns to its starting point. By the estimate on the 
expected hitting time (3.2) of [23], it follows from (|6.9p that 

hm J-Q^^lHixf")] = -, and, analogously, for j ^ 1, (6.10) 
1™ = - and hm l^Q^jHiA^' ^ ^ 



because, by Lemma [Ol the escape probability Vd equals the capacity between the 
origin and infinity in Z''. Note that A(0) = 1 and that the capacity is computed 
with respect to the counting measure. 

Define 5* to be the stopping time given by the first visit to after the first visit 
to Xi , after visiting A^^j^. By the strong Markov property, Q^[S'] is equal to 



Hence, by (lOO)) . 



We are now in a position to prove the lemma. Fix an arbitrary sequence in 
BiA^jJ^y. Since (Q^„ [#{visits to before H{A^j,)}] = Q^„[#{visits to < 
before H{A^j ^)} l{H{x^) < i/(A]^^ J}], by the strong Markov property, 



[#{visits to x^ before H{Afj^^)}] 
J^„[#{visits to x^ before H{A^, ,)}] 



}^^[S] = [#{visits to xf before 5}] 



To estimate the numerator, observe that, by [Ij, Chapter 2, Proposition 3 and 
Lemma 7, 

Nd 

= Qf«[#{vis. toxf bef. i?«f,i)}] + Q^«[#{vis. toxf bef. F(y^)}] 

= Qf„[#{vis. toxf bef. i?«i)}] + ^{Q^«[i/(?/^)]+Qf«[i?(xf)]}. 
Hence, the right hand side of (|6.12p can be written as 

A^-'' [^] - Qf« [i? (y~)] - Qf« [i? «)]) 

Q^„[#{visits to xf before i/«) o e^^^.^^^ + i?«,,i)}] ' 

Let 5" be the stopping time H{A^.j -^) + H{x^) o Oj^f^^N By Lemma 7, Chapter 2 
in [T], the denominator is equal to 7V~''Q^jv[S"]. Hence, by the strong Markov 
property on H{A^j the previous ratio is equal to 







-Qf«[i?(y~)]-C 










= xf]Q^AHix^)]) 



QUENCHED SCALING LIMITS OF TRAP MODELS 



27 



By (|6.10p . (16. lip , this expression converges to M ^ as N 1 oo. Since the sequence 
G B{A^ , InY is arbitrary, we are done. □ 

For a subset F of T^, let f{F) be the return time to F for the discrete time 
random walk {Yj^ : k > 0}: 

fiF) = mi{k > 1 : 1/ e F} . 

For subsets A, B of such that A n i? = 0, let capyjv (A, B) be the capacity 
between A and B induced by the process : 

capy„(A,i?) = infl E[/(2^)-/(^)]'' 

where the infimum is carried over all functions / : M such that f{x) = 1 for 

all X in A, f{x) = for all x in B. 

Corollary 6.7. Under the same conditions of Lemma \6.6\ for j ^ 1, 

hm Q^^[i/(xf)<f«,)] = ^. 

Proof. Since 

lim Q^4H{dB{xf,lN))<T{^f)] = Vd, 
the result follows from the strong Markov property and Lemma 16.61 □ 

Corollary 6.8. If In 1 oo and x^ £ Tf^, are such that d{x^ ,y^) > In, 

lim capy«(a;^,y^) = ^- (6.13) 

Proof. The corollary is a direct application of Lemma 16.21 and Corollary 16.71 □ 

6.3. Random walk on for d ^ 2. We present similar results to the ones stated 
in the previous subsection, but which also hold in dimension 2. Here, however, we 
need to impose that the distances In, between the very deep traps and the starting 
point of the walk, grow close to linearly. 

Although Lemma 16.91 below also holds for d ^ 3 and could be used in place of 
Lemma |6^ we keep both results since they are based on different arguments. Let 
In be an increasing sequence such that In/N and In/N'^ oo for every a <1. 
In this way, 

lini k = , lim iHli^ = 1 . (6.14) 

N^aa N N^oo log N 

In this subsection it will be more convenient to work with the distance d2{x,y) 
in given by N times the Euclidean distance between x and y in T"^. 

Lemma 6.9. Consider a sequence of sets A^j ~ {x^ ,...,x^[} C such that 
d2{xi,Xj) ^ In (0 ^ i < j ^ M) for some sequence {In '. N > 1} satisfying (I6.14p . 
Then, 

e[H(xf)<i/«jl-4 = 



lim sup 



M 
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Since the result only concerns the first point in A^j to be visited, we can suppose 
that the process {Yj. : /c > 0} is a lazy random walk in T^, i.e. with probability 
one half Y does not move, otherwise it jumps uniformly to one of its neighbors. 
Before going into the proof of Lemma [5?^ we collect some properties of the hitting 
and mixing times of Y . 

Recall that || • || denotes the total variation distance between two probability 
measures. The following bound on the mixing time on the torus follows, for instance, 
from Corollary 5.3 and equation (5.9) in [T7] . 

hm limsup||Qf[r^Ar2 =.]-7r^(.)|| =0. (6.15) 

Of course, the same result holds for any sequence {tAr : > 1} which increases to 
oo faster than iV^. 

We claim that for every /? > and x'^ ,y'^ £ such that d2{x^ ,y^) ^ In 

lim Q^^[Hix^)^PN^] = 0. (6.16) 

Since for d ^ 3 this statement follows from (j6.6|) . we concentrate in the case d = 2. 
Recall that Y stands for the simple random walk on I? (with law P), and denote 
by (j>N ■ 2^ ^ TT^ the canonical projection. In view of the invariance principle in 
Z^, it is enough to prove that for every i? ^ 0, 

lim FyN visits x^ before Y exits B{y^ , RN)] = , 

A^ — *oo 

where £ is such that (fiNiy^) = G T^- This follows, for instance, from 
[7], (225) since for every point x'^ g B{y'^ , RN) such that (j){x^) = x^ , we have 
\y^ — ^ In and, moreover, the number of possible choices for x^ is bounded 
uniformly on iV ^ 1. This establishes (|6.16p for d — 2. 
A straightforward consequence of (|6.16p is that 

lim Q^r^iHix^) ^ (3N^] = (6.17) 

for any sequence {x^ : iV > 1} in T^. 

Now we can state the hitting time estimate we will use during the proof. Consider 
the scales 

^N^logN lid ^2, 
N"^ if d ^ 3 . 



We claim that for every d ^ 2, 

lim lim sup sup Qy[H{x) < jh%] = . (6.18) 

Indeed, fix 5 > 0. By (|6.16p . for any /3 > 0, for iV sufficiently large, and for any 
X, y such that d2{x, y) > In, 

[H{x) < < 6 + Q^[H{x) > PN^ , H{x) < jh%] . 

On the set H{x) > /3N'^, H{x) = fiN"^ + H{x) o dpN^. Therefore, by the Markov 
property at PN'^, the second term is less than or equal to 



[Hix) < ^hl,] 



By (|6.15p . for f3 large enough, this expression is bounded by 

S + Qf« [H{x) < jh%] . 
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The result follows now from Theorem 2.1 of [23] and Lemma 16^ 

Finally, we claim that for every positive 7, the probability to hit a very deep 
trap before 7/1^ is bounded away from zero in N. More precisely, for every 7 > 0, 



N—KyD 



< 1. 

The claim above also follows from Theorem 2.1 of 123 and Lemma 



(6.19) 



Proof of Lemma \6.9l Our strategy is to consider several consecutive attempts to 
hit one of the points in A^,j. For 7 > and a positive integer L, define the times 



a, = i{LN^+jh%) 



i{LN^+jh%) + LN^ 



for i > 0. Intuitively, for each i, we use the intervals [a^, bi] to approach equilibrium 
measure tt^ and the intervals [6.^,024.1] to attempt to hit the set A^j. 
Let 



Rl n ■— max 



Sl.n := sup 



By symmetry, the maximum is irrelevant in the definition of Rl.n, and, by (|6.15p . 
Rl,n vanishes as iV | 00 and then L "f 00. By (j6.16p . (j6.17p . Sl.n an Tl^n vanish 
a.s N 1 00 for every L > 1. 

For < s < t, define the random variable Js,t, which takes the value if the 
set A^j = {^^fi • • ■ i^mI is ^ot visited between the times s and t and otherwise 
Js^t = 1,...,M according to the index of the first point in A^j visited in this 
interval. 

The proof of the lemma is divided in two parts. We first claim that for every 
7>0, 



lim sup sup 



N 



[Jbo.ai - 1] 



t,N 



. (6.20) 



Note that this expression does not depend on L, but only on 7 and N . Then, we 
prove that 

M 



lim lim sup 



= . 



(6.21) 



Clearly, Lemma [Ql follows from and (p?^ . 

The proof of (|6.20p relies on three estimates. Consider the event 'Dp — [Jat.bi 7^ 
for some i = 0, . . . , F — 1], F > 1. This event indicates that some very deep trap 
is visited in one of the "mixing" intervals [aj,bj]. We claim that 



sup 



'y [Sf] < FTl.n + Sl.n + Rl.n 



(6.22) 



Fix a site y in B{A^_j,lj\!y and decompose the event Dp according to whether 
the set A^j has been attained before time LiV^ or not to get that 



tiA 



[Dp , H{AZ) > LN'] 
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The first term is bounded by Sl.n- On the set iJ(Aj^) > LN'^, the event S);- is 
equal to Ui<i<_F_i{ Ja^.b^ 7^ 0}. Hence, by the Markov property at time bo, the 
second term on the right hand side of the previous inequality is bounded by 



En 



F-2 



< Rl. 



N 



i=0 



Since tt is the stationary state, the second term is smaller than or equal to {F — 1) 



Ar[-ff(A^) < LN ]. In conclusion, we proved that 



1^' 



N 



which is exactly (|6.22p . 

Let > 1, be the event that no site in A^j has been visited in the "hitting" 

intervals a^+i], < i < F - 1: = [Jb,,ai+i = for aU i = 0, . . . , F - 1]. We 
claim that 



sup 



[€f] - =0]^ 



< FRl. 



N 



(6.23) 



Fix a site y in B{Aj(j, InY- By the Markov property, 

F-2 

Q^I^f] = [ n { Jft.,a. + , - 0} ^ [Jt,,a, = 0] 

4=0 " ^ 

Applying the Markov property at time 60 — LN'^, this expression can be written as 

<[ n iJb„a,^, = 0}](Q^„[Jo,ai-bo = O] + Rn 
i=0 

where the remainder Rn is absolutely bounded by Rl.n- Since tt^ is the stationary 
state, Q^N[Jn.ai-bo = 0] = Q^N[Jbo,ai =0]. We proceed by induction to derive 

Let Sjp, F>l,he the event {H{xi) < n n We claim that 



sup 



[Jbo.a, ^ 0] 



< F^Rl,N + Q^N[Jbo,a,=0] + Q^I^f] 



(6.24) 



Clearly, 



F-l j-l 

= S)^ n y { fl {Jb.^a,^, = 0} n {Jb^^a,^, - i}} . 



j=0 k=0 

Repeating the arguments which leaded to (|6.23p . we get that for each < j < F — 1 



and any site y in B{A^,j, InY, 



k=0 



[Jbo,ai — 0]''Q^n \Jbo,ax — l] + R-N , 
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where the remamder Rn is absolutely bounded by [j + 1)Rl,n- The value of the 
remainder Rn may change from line to line below. Summing over j, we get that 

F-l 

^v[^f] = Qfiv [Jfc„,ai = 1] E C« [\>■a^ = O] ^' + Rn , 

i=o 

where the remainder Rn is now absolutely bounded by F'^Rl,n + Qyl'^p]- This 
expression can be written as 

=1] , 

Kn 



l-Qf«K.ai -0] 

for a remainder Rn absolutely bounded by F'^Rl^n + Q^n [Jbg.ai = 0]^ + Q^[X)f], 
which is precisely (|6.24p . 

By the estimates (I02)) . ((03t . (|6T24| . the supremum in (|6?20| is bounded by 

2FTl,n + 2Sl,n + ^F^Rl^n + 2Qf« = O]^ 

for every F > 1, L > 1, 7 > 0. By (pT^ . Qf„[Jfco,ai = 0], which does not depend 
on L, is bounded above by a constant strictly smaller than one. Hence, as iV | 00, 
and then L ] 00, Rl,n, Sl.n and T^jv vanish. It remains to let F t 00 to conclude 
the proof of ((OOl) . 

It remains to prove (|6.2ip . Decompose the event {Jba,ai = 1} according to the 
event that at least two sites in have been visited in the time interval [bo, ai] to 
get that 

iQ^AJbo.a, = 1] - Qf«[Y[fc„,„,] n - {xi}] I < Qf« [#(r[b„.a,] n <) > i] . 

In this formula, yfhg.ai] stands for the sites visited by the random walk in the interval 
[5o,ai], yjbo.ai] — {Yt ■ bo < t < ai}, and ^A for the cardinality of A. By similar 
reasons, 

\QnAy[bo..a,] n AZ = {Xi}] - Qf„[xi e Y[fc„,,,]] I < Qf„ [#(y[bo,a,] H > l] . 

Therefore, 

iQ^AJbo.a, = 1] - Qf«[xi e Y[fc„,„,]] I < 2Q^„ [#(r[6o,a,] n <) > i] . 

Note that the probability Q^iv[y € ^[60, ai]] does not depend on y by symmetry. In 
particular, it also follows from the previous arguments that 

\Q^4Jbo,a, ^ 0] - MQ^4x, e y[,„,,,]] | < (m + i)Qf„ [#(F[bo.a,] n ^j^,) > i] 

so that 

\QUJbo..a^ = 1] - (l/M)(Qf„[Jfc„,„, ^ 0] I < 4Q^„ [mYM n O > 1] . 
Equivalently, 



Jf« [jbo..a, = 1] 1 



[Jbo.ai 7^ 0] M - Qf« [j„o,ai 7^ 0] 

By the strong Markov property, the right hand side is bounded above by 

max Q^Ji7«,,) < . 

By ()6.18p . this expression vanishes as iV t 00 and then 7 i 0. This concludes the 
proof of the lemma. □ 
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We also need to estimate in dimension 2 the probability that the random walk 
escapes from a deep trap. More precisely, 

Lemma 6.10. Assume that d ~ 2 and let {Ij^ : N > 1} be a sequence satisfying 
(|6.14p . Then, for any sequence {y^ G : N > 1}, 

lim log7V(Qf«[i/(S(y^,M^)<f(y~)] = J- 

Proof. The number of visits that the random walk performs to y^ before exit- 
ing B{y^ ,1m) is a geometric random variable, with failure probability given by 
Qyjv[H{B{y^ < T'iy'^)]- The inverse of this probability is equal to the ex- 
pected number of visits to y^ before exiting B{y^ ,1^)- By [1^, Theorem 1.6.6, 
the expected number of visits, denoted by GB{yi^ jj^){y^ ,y^) in the notation of 
Green's functions, is given by 

GBiy-M)iy''^y'') = ^log^ + K + o{N-') , 

for some constant K in R. The result follows from this estimate. □ 



Corollary 6.11. Assume that d — 2. Under the hypotheses of Lemma \6.9\. for 

Jim^logiVQ^„[i7(xf)<f«,J] = 

Proof. This result follows from the strong Markov property and Lemmas 16.91 and 
QUI □ 



Corollary 6.12. Assume that d = 2 and let {Rn '■ N > 1} be a sequence such that 



Rn 1 oo. Let ,y^ £ T^, such that d{x^ ,y^) > Rjy. Then, 



lim logiVcapy«(a;^,?/^) = ^- (6.25) 

AT— *oo 4 

Proof. The corollary is a direct application of Lemma 16.21 and Corollary 16.111 □ 

6.4. Metastability of the trap model in dimension d > 3. Recall that we 
denoted by Vd the probability that a nearest-neighbor, symmetric random walk on 
Z'^ never returns to its starting point. As in the previous subsections, we denote 
by Yj^ the discrete time random walk on the torus T^, inducing the law on 
i?(Z_|_, T^). The proof of Theorem l2.3l is divided in two parts. In Proposition l6.13l 
below we show that the trace process converges and in Corollary 16 . 151 that the time 
spent outside A^ is negligible. 

Proposition 6.13. Fix M > 1 and T > 0. As N ^ oo, the process {Xf '^^ : 
< t < T} converges in distribution to the Markov process on {1,...,A/} with 
generator Lm given by 

M 

Proof. Fix M > 1 and denote by r^.M ■ {1, ■ • ■ , M} x {1, . . . , M} M+ the jump 
rates of the trace process {X^''^'^ ■ < t}. By (|6.1|) . for j ^ i, 
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where again A^jj — {xi , ■ ■ ■ , Xj^_i,Xj'^i,x^j}, I < j < M < iV*. To compute this 
probabiHty, we need only to examine the discrete skeleton Markov chain: 



N 



Since xj^ converges, as TV t oo, to Xj, 1 < j < M, and since mini<i^j<M ll^^i — > 



0, by Corollary EJl 

Vd 

TV— foo 



M 



r Hixf)<riAZJ_ 

Hence, for j ^ i, riy,M{hj) converges, as t oo, to {vd/Mwi), because 
converges to lii. This concludes the proof of the proposition. □ 

To examine the time spent by the random walk {JVf : < < < T} on \ A'j^j, 
denote by cap^ the capacity associated to the process . Of course, for any two 
disjoint subsets A, B of T^, 

capjv(A,B) = -^;^pi;^cap^«(A,B) . (6.26) 

For X ^ y, in Tff, denote by P^'^ the probability measure on the path space 
D{R+,T% \ {x}) induced by the trace of {X^ : t > 0} on \ {x} starting from 
y. Expectation with respect to P^^^ is denoted by E^'^. 

Lemma 6.14. We have that 



lim limsup max max AfE„ ^ HlA^f A 

Af^OO JV^oo 1<J<M y:\y-x!^\ = l 



. 



Proof. Fix 1 < j < M and y ~ xf . By Lemma 16.11 the expectation appearing in 



the statement of the lemma is equal to 



1 



A- ) ^ -''^)^'^\ii^y) < H{AZ,)] ■ (6.27) 



z=ixf 



By (|6.26p . the denominator is equal to 



■capyiv(y,A]^^^^) > — ^ capyiv (y, ) 



In view of CoroUarv 16.81 this latter expression is bounded below, uniformly in N , 
by a strictly positive constant. 

To estimate the numerator in (j6.27p . we need only to examine the discrete skele- 
ton Markov chain because Pf [if (y) < H{Al^ -)\ = [ff(2/) < li{A%^^)\. Fix a se- 
quence {In --N > 1} such that 1 « In « N and let Bn = {z : d{z, A^j) < 
(n} \ AZ, CN = {zeT%: d{z, AZ) > In}- Since Q^[H{y) < i?«,,,)] vanishes 
on the set A^ j , 

^ .^(z) [H{y) < i/«,)] = [^(y) < ^«..)] 

+ 5] ^^(z)Qf[ff(y)<i/(A]r,,,.)] • 

The first term on the right hand side is bounded by ^^{Bn), which vanishes as 
N ^ oo because £n « N. Since £n ^ 1, by Lemma WM as N ^ oo, the second 
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term converges to M ^[1 — VK(T'^) ^ J2i<j<M ''^j]- The expression inside brackets 
vanishes as M "f oo by definition of the sequence {wi}. This proves the lemma. □ 

Recall that An,m = ^%\ ^m- 

Corollary 6.15. For every t > 0, 

lim limsup max E^nItJ^"''^'] = 0. 

M^oo N^oc l<j<M L c J 

Proof. Fix M > 1 and 1 < j < M. Consider the stochastic process Z^'^^ with 
state space defined as 

where a{t) := sup{s < t : e ^m}- Hence, during an excursion in An,m by 
X'^ , the process ' stays at the last visited site in A^^. 

For a path uj G D{M.+ ,Tff) performing infinitely many jumps, denote by T„(a;), 
n > 0, the jumping times of u: to(w) = and 

Tn{uj) := inf{t > T„_i(cj) : uj{t) ^ cj(T„_i(a;))} . 

Let 

T„{uj) := T„{uj) - Tn-i{uj) , n>l, 
and let Nt be the number of jumps up to time t: 

Nt{Lu) := sup{j > : t,(w) < t} . 

The process X^''^^ , defined on the path space _D(R+,y4j^), can be thought as 
a process on Z?(R+,T^). Couple the processes X^''^^ , Z^'^^ forcing them to visit 
the same sequence of sites. By Lemma [6.141 and the proof of Lemma 4.4 in [2 , for 
every K > 1, 

lim limsup max [ r^Af (^^'*'^) - tkm(^^'^'^) ] = 0. (6.28) 

Af^oo ^^j^^i J 

Set Nt Nt{Z^'^^), f„ and T„ T„(l^'*^). Fix 1 < j < M. 

Under P^„, 



Nt + l KM 

Tt^"-"' < i A ^ (f„ - r„) < l{Nt > KM} t + Y^iTn- 



n—1 n—1 

for any positive integer K. Therefore, 

E^ivK^"'"'] < tP^4Nt>KM] + E^4TKM{Z''^^')-rKM{X''^^')] - 
j j j 

By (I6.28p . the second term vanishes as t and then M f oo. It remains to 

prove that 

lim limsuplimsup max P'^j^[NtiZ^-^^) > KM] = 0. (6.29) 

K^co M^oo N^oc i<j<M 

Since NtiZ^'^^) < Nt{X^'^'), - a.s., 

P^„[A^t(Z^'*0 > KM] < Pf«[iVt(X^^*0 > KM] . 

j j 

Fix M >1 and 1 <j < M such that 

limsup max F^^[Nt{X^'^'') > KM] = lim sup Pf« [iVt(l^'*^) > i^M] . 
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Since [Nt > KM] is a closed set for the Skorohod topology on D{M.+ , A^,j), since 
Nt{X'^''^'^) has the same distribution as Nt{X^''^'^) and since, by Proposition [OSl 
j^N,M converges in distribution to Z^'^ , 

hm sup [Nt {X^'^^) > KM] < max Pk [Nt (Z^) > KM] , 



l<k<M 



where Pk is the distribution of the process Z starting from k. 

To estimate the right hand side, we compare Nt{Z^'^) with a counting process Ct 
in which we replace the holding times T„ by if Z*^(r„_i) 1. In other words, let 
Go := Co be the number of times the process jumped before hitting 1 for the 
first time. Since jumps from any site uniformly to all others. Go is a random 
variables with geometric distribution: P[Go = n] = {1/M)[{M - l)/M]"-\ n > 1. 
When hitting 1, as Z^'^ , the process Ct stays there for a mean wi/vd exponential 
time. At the end of this exponential time, Ct jumps from Go to Go + Gi, where 
Gi stands for the number of jumps performed by Z^^ before hitting 1 again. 

By construction, Nt{Z^'^) < Ct for all t > and Ct = J2o<j<Nt ^i' where 
{Gj : J > 0} arc i.i.d. random variables with geometric distribution: P[Gi = 
n] = (1/M)[(M — 1)/M]"'~^, n > 1; and Nt is a Poisson process with rate Vd/wi, 
independent of the sequence {Gi}. In particular, E'p^^ [Nt{Z^^)] < M[l + (twi/vd)]- 
This proves (|6.29p and the corollary. □ 

Proof of Theorem \2.3[ Theorem 12.31 follows from Proposition 16.131 and Corollary 
EH □ 

Proof of Theorem \2.5\ Denote by {Zt : i > 0} the A'-process with parameters 
{wi/vd : j > 1} and c — Q. Using independent exponential and Poisson random 
variables, we may define in the same probability space (f2, JF, P) processes {X^'^^ : 
i > 0}, {Zf^ : t > 0} and {It ■ t > 0} which have the same distribution as 
{Xf '^^ : i > 0}, {Zf' : t > 0} and {Zt : t > 0}, respectively. Fix a common 
starting point j for all processes and T > 0. By [T^l Lemma 3.11], {Zj^ : < t < T} 
converges a.s., as M t oo, to {Z^ : i > 0} in the Skorohod metric. On the other 
hand, by Proposition 16. 131 if ds stands for the Skorohod metric on D([0, T], N), for 
every M >1 and e > 0, 

lim Pfe(X^'*^Z*0 > el = . 
In particular, there exists a strictly increasing sequence {Nlj : M > 1}, such that 

for all TV > N'^.j. Hence, by the triangular inequality, for any sequence {Nm : M > 
1} such that Nm > Nlj, 

lim P^ds(X""^*^Z) > e] = . 

for any e > 0. The sequence {l*^ : TV > 1}, defined as the inverse of {Nlj : M > 1}, 
fulfills the requirements of the first part of Theorem 12.51 

To prove the second statement of Theorem 12.51 fix M > 2 and observe that 
^ HiA'^if) + Tt^"-"' o e{H{AZ)) provided iN > M. In this formula. 
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0{s) : D{R+,T%) D{R+,T%), s > 0, stands for the time shift by s of a path oj: 
{e{s)uj){t) =Lo{s + t),t>0. Therefore, 

max E^^fr^"''"] < max EWH(AZ)] + max fT,^"'"] 

provided M < i^. By Corollary I6.15[ the second expression converges to as 
iV t oo, M t oo. By definition of H{A^), the first expectation on the right hand 
side vanishes for 1 < j < M. For M <j < l^, by ([13]) with F = E = T%, 

^ cap(xf ■ 

The denominator is bounded below by min^g^jd^ maxi<j,<jv/ cap(z, a;^). Since \\X\ — 
X2\\ > 0, by Corollarv l6.81 this expression is bounded below by a positive constant, 
uniformly in M and N. This proves the second statement of Theorem 12.51 because 
ly'^ [An^m] vanishes as iV "f oo, M "f oo. □ 



6.5. Metastability for d = 2. In this subsection, we adapt to dimension 2 the 
results presented in the previous subsection. Most of the proofs are similar to the 
case c? > 3. 

The main difference with respect to dimension 3 is that the process is speeded 
up by logiV. Recall from Section [5] that we denote by {Xf : t > 0} the random 
walk with generator Cn, defined in (|2.3p . speeded up by log TV. Moreover, P^, P^, 
X £ T^, stand for the probability measure on D(R-|_, T^) induced by the processes 
{Xf : i > 0}, {Xf : t > 0} starting from x. 

Proposition 6.16. For a fixed M > 1, and T > 0, the process {xf : t T} 
converges in distribution, as N '\ oo, to the Markov process in {1, . . . , M} given by 
the following generator 

Proof. As in the proof of Proposition 16.131 we use (|6.ip to write the jump rates of 
Xf^'*^ in terms of the excursion probabilities between the very deep traps: 



Note the factor log N which appears because the generator C m is multiplied by this 
constant. 

Consider a sequence {In '■ N > 1} satisfying (|6.14p . Use the strong Markov 
property on H{B{xf ,ImY) to obtain 



)] 



l{i7(i?(xf,Z^r)<f(xf)} 



[H{xf) < H{A 



M,j] 
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Therefore, 



logiVQ^^[H(xf)<f(A]v,)]-:^ 



< 



2M 

„N 



H{B{x^Mr)<T{x^)\ - 



sup 

z&dB{xf ,l„) 



1 

M 



By Lemmas 16.91 and 16.101 these expressions vanish as N 1 oo. Since W^n converges 
towards Wi, 1 < i < AI, we are done. □ 



Recall the definition of the measures 
previous subsection. It corresponds to the trace on \ {x} of the process {X/ 
t >0}, which has not been speeded up. 



'y''^, X ^ y e T^, introduced in the 



Lemma 6.17. In dimension 2, 



lim lim sup max max 

M^oo N^oc ^<3<M y:\y-xf \ = l logiV 







Proof. The proof of this result follows the same argument as in Lemma 16.141 One 
only notes that the denominator of (|6.27p is now multiplied by log N which allows 
us to use Corollary 16.121 in place of Corollary 16.81 The argument to bound the 
numerator is also the same. However, one should choose a sequence {^n '■ N >1} 
satisfying (|6.14p in order to apply Lemma 16.91 □ 



For X 7^ ?/, in T^, denote by P^^^ the probability measure on the path space 
D(K+,T5^ \ {x}) induced by the trace of {Xf : i > 0} on \ {a;} starting from 
y. Expectation with respect to P^'^ is denoted by E^'^. The difference between 
P^'"^ and P^'"^ is that the first probability measure is associated to the random 
walk speeded up by logiV. Therefore, for every subset A of \ {x^}^ 



N.: 



[H{A)\ = 



log TV 



[H{A)\ 



In particular, it follows from the previous lemma that 
lim lim sup max max ME„' ^ \H{A 

M^oo JV^OO 1<]<M y:\y-xf \=1 

Corollary 6.18. In dimension 2, for every t>Q, 

lim lim sup max E^'^w fe'^"'''] 

Af^oo jv^oc- l<j<A/ '^i 



M,j) 



(6.30) 







Proof. The argument is identical to the one in d ^ 3 presented in CoroUarv 16.151 
We just use (|6.30p and Proposition 16.161 instead of Lemma 16.141 and Proposition 
16.131 At the end of the proof, the rate of the process 7Vt(Z*^) is replaced by tt/ {2wi), 
but its exact value is superfluous. □ 



Proof of Theorem \2.4\ The proof is a direct consequence of Proposition 16.161 and 
Corollary iH □ 
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6.6. Dimension 2 with no acceleration. We prove in this subsection that in 
dimension 2 the trap model with generator ()2.3|) starting from a very deep trap 
does not move. Hence, on the order 1 scale, the random walk does not move and 
on the scale log iV it converges to the if-process in which all the geometry is wiped 
out. 

Proposition 6.19. For every j > I, every t > and every sequence {£n : .^V > 1} 
such that £n T oo, 

lim ¥^.[\X^it)~xf\>£N] = 0. 

Proof. Fix j > I and a sequence {£n ■ N > 1} such that £n oo. Following [5], 
denote by Sn : Z+ ^ M the clock process: S'Ar(O) = 0, 

fc-i 

SN{k) - ^e,M^^(r^(z)) , k>l, 

2=0 

where {Y^{i) : i > 0} is a nearest-neighbor, symmetric, discrete time random 
walk on starting from ; {ti : i > 0} is a sequence of i.i.d. mean one, 
exponential random variables, independent from the Markov chain {Y^{i)}-, and 
W^{x) = W^, X eT%. Denote by Tat : K+ ^ K+ the inverse of S_ 



N- 

TN{t) = supjfc : SN{k) < t} . 

Clearly {X^(t) : t > 0} has the same distribution as {Y^{TNit)) : t > 0}. 
Hence, 

F^^[\x''{t)~xf\>£^] = P''[\Y''{T^{t))-xf\>£^] 

< max \Y^{k)-xf\>£j,] + P^[Tj,it)>rj,] , ^^'^^^ 

for a sequence rjv such that 1 « r^ « £%■ 

We estimate separately the expressions on the right hand side of (|6.3ip . Since 
Tjv is the inverse of Sn, {TN{t) > r^} = {S'jv(''jv) < i}- In particular, 

rjv-l 

P^[r^(i) > tn] < P^K- e. = xf} < t] , 

i=0 

because SN{k) > W^, Eo<^<k-l ^^ HY^{i) = xf}, W^, > w,. Since {r^(fc) : 

3 — — j 

fc > 0} starts from x^ and the two-dimensional random walk is recurrent, the 
previous probability vanishes as N ] oo because rjv T oo. 

On the other hand, since Y^ [k) — x^ is a bi-dimensional martingale, by Doob's 
inequality, 

max |r^(fc)-xf|>M < 
which vanishes as N ^ oo. This concludes the proof of the proposition. □ 
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